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ABSTRACT 

Detection of gamma rays from dark matter annihilation in the galactic center is 
^ ' one of the feasible techniques to search for dark matter. We evaluate the gamma 

ray flux in the case that the dark matter has an electroweak SU{2)l charge. Such 
^ • dark matter is realized in the minimal supersymmetric standard model (MSSM) 



when the lightest SUSY particle is the Higgsino- or Wino-like neutralino. When 
the dark matter is heavy compared to the weak gauge bosons, the leading-order 
calculation of the annihilation cross sections in perturbation breaks down due to a 
threshold singularity. We take into account non-perturbative effects by using the 
non-relativistic effective theory for the two-body states of the dark matter and its 
SU{2)l partner (s), and evaluate precise cross sections relevant to the gamma ray 
fluxes. We find that the annihilation cross sections may be enhanced by several 
orders of magnitude due to resonances when the dark matter mass is larger than 
1 TeV. Furthermore, the annihilation cross sections in the MSSM may be changed 
by factors even when the mass is about 500 GeV. We also discuss sensitivities to 
gamma ray signals from the galactic center in the GLAST satellite detector and the 
large Air Cerenkov Telescope arrays. 
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I Introduction 

Recent cosmological observations determine precisely the mean densities of mat- 
ter and baryon in the Universe and existence of non-baryonic dark matter is 
estabhshed. Weakly interacting massive particles (WIMPs) are considered to be 
good candidates of the dark matter 12.j. They act as the cold dark matter in the 
structure formation of the universe. High resolution A^-body simulations show that 
the cold dark matter hypothesis explains well the structure larger than about 1 Mpc 
jSj . On the other hand, fundamental problems, such as (i) the constituent of the dark 
matter and the origin in the thermal history and (ii) the dark matter distribution in 
the galactic scale, are not still resolved. It is important to detect the dark matter in 
direct or indirect methods in order to answer the questions. 

Many detection methods have been proposed and some of the experiments are 
now operating. Among those, the detections of exotic cosmic ray fluxes, such as 
positrons, anti-protons and gamma rays, are feasible techniques to search for the dark 
matter particles In particular, an excess of monochromatic (line) gamma rays 

due to the pair annihilation would be a robust signal if observed, because diffused 
gamma ray background must have a continuous energy spectrum 0. The GLAST 
satellite detector [TU| and the large Atmospheric Cerenkov Telescope (ACT) arrays, 
such as CANGAROO III [TT], HESS [T^, MAGIC ^ and VERITAS dH, can search 
for the exotic gamma rays from the galactic center, the galactic halo, and even from 
extra galaxies. 

In this paper, we discuss SU(2)l non-singlet WIMPs and the gamma ray fluxes 
from the galactic center due to the pair annihilation. We refer to such dark mat- 
ter as electroweak-interacting massive particle (EWIMP) dark matter ^21- Due to 
their SU(2)i non-singlet nature, EWIMPs have interactions with the SU(2)i gauge 
bosons such as W and Z bosons. If EWIMPs have a vector coupling to Z boson, 
the current bound obtained from direct dark matter searches through their spin- 
independent interaction is stringent [16^. This means that the EWIMP dark matter 
should be a Majorana fermion or a real scalar if the mass is around 1 TeV. In this 
paper, we consider the former case. We especially study triplet and doublet EWIMP 
dark matters, which are neutral components of an SU(2)i-triplet fermion whose hy- 
percharge is zero and of a pair of SU(2)j;^-doublet fermions with hypercharges ±1/2, 
respectively. When the EWIMP mass is around 1 TeV, the thermal relic abundance 
may be consistent with the cosmological observation. 
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The EWIMP is realized in the minimal supersymmetric standard model (MSSM) 
when the lightest SUSY particle (LSP) is the Higgsino- or Wino-like neutralino [T7] . 
Wino is the superpartner of the SU(2)i gauge boson, and Higgsino is that of SU(2)i 
doublet Higgs bosons. The thermal relic density is too low if the LSP mass is 
smaller than 1 TeV. However, decays of gravitino or other quasi-stable particles may 
produce the LSPs non-thermally so that the relic abundance is consistent with the 
cosmological observation. While the LSP with the mass about 1 TeV may lead to 
the naturalness problem, such possibilities are discussed in the split SUSY scenario 

The line and continuum gamma ray fluxes from the EWIMP dark matter anni- 
hilation in the galactic center are proportional to the cross section to two photons 
and those to other modes whose final states fragment into vr^s, respectively. The 
leading-order cross sections in perturbation have been calculated by many authors 
|19j . However, if the EWIMP mass is large compared to the weak gauge boson 
masses, usual perturbative expansion for the annihilation cross sections cannot be 
applied ^^l- This can be seen in the violation of unitarity of the one-loop annihi- 
lation cross section into two photons. The fact comes from the degeneracy of the 
EWIMP and its S\J{2)l partner(s) in mass and the non-relativistic motion of the 
dark matter in the current universe. The transition between EWIMP and its partner 
pairs is induced by the t-channel weak gauge boson exchange. When the EWIMP 
mass is much larger than the weak gauge boson masses, the weak interaction behaves 
as a long-range force. The wave functions of EWIMP and its partner pairs are mod- 
ified from plane waves at the non-relativistic limit, and the mixing between those 
states is enhanced. This phenomenon is related to so-called a threshold singularity, 
and we have to consider the effects of the long-range force on the annihilation cross 
sections for reliable calculation. 

In this paper we work in the non-relativistic effective theory for EWIMPs. Non- 
relativistic effective theories are often used in calculations of the threshold pro- 
ductions of heavy particles, the quarkonium mass spectrums and so on |2I]. In 
this technique, we can factorize short-distance physics, such as pair annihilation, 
from long-range effects on the wave functions due to the optical theorem The 
long-range effects are evaluated by solving the wave functions under the potential. 

We found that the annihilation cross sections may be enhanced by several orders 
of magnitude compared to the leading-order calculation in perturbation when the 
EWIMP mass is larger than about 1 TeV. The mixture of the pairs of EWIMPs and 
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the SU(2) L partners can form a bound state whose binding energy is close to zero, 
and it contributes to the annihilation cross sections in the non-relativistic limit. The 
enhancement of the cross sections originates from the resonance by the bound state. 
Furthermore, the annihilation cross section to two photons, which is suppressed by 
a loop factor in perturbation, becomes comparable to those to the other modes 
around the resonance. As a result, the continuum and line gamma ray fluxes from 
the galactic center due to the EWIMP annihilation are enhanced. The indirect dark 
matter searches by the large ACT detectors, which have sensitivities to TeV-scale 
gamma rays, may be promising, if dark matter is a TeV EWIMP. 

It is also found that the non-perturbative corrections to the cross sections are 
sizable for the triplet (doublet) EWIMP even when the EWIMP mass is about 500 
(1500) GeV. Thus, the correction should be taken into account in the evaluation of 
the gamma ray fluxes in the MSSM, especially when the LSP is Wino-like. 

This paper is organized as follows. Wc flrst summarize the properties of EWIMPs 
and discuss the threshold singularity in the EWIMP pair annihilation in the next 
section. In Section III the non-relativistic effective actions for the triplet and the 
doublet EWIMP pairs are derived. In Section IV the cross section formula is obtained 
using the optical theorem and the effective actions. While the obtained annihilation 
cross section to two photons in our formula is reduced to the one-loop result in the 
perturbative expansion, it also satisfies the unitarity bound in the limit of an infinite 
EWIMP mass. The one-loop cross section does not satisfies the bound. 

In Section V some numerical results for the annihilation cross sections are pre- 
sented. The fitting formulae for the annihilation cross sections are also derived from 
the numerical results. The resonance behaviors of the cross sections are studied us- 
ing a toy model in which the electroweak Yukawa potentials are approximated by a 
well potential. In Section VI we evaluate the gamma ray fiuxes from the EWIMP 
annihilation in the galactic center and discuss the sensitivities in the future exper- 
iments. In Sections V and VI the cross sections and the gamma ray fluxes for the 
Wino- and the Higgsino-like ncutralinos are also evaluated in the wide range of the 
MSSM parameters. Section VII is devoted to summary of the paper. 

II Properties of EWIMPs and Threshold Singularity 

In this section the mass spectrums of EWIMPs and the SU(2)i partners and their 
low-energy interactions are summarized. We also discuss the threshold singularity 
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in the non-relativistic EWIMP pair annihilation. The perturbative expansion of the 
annihilation cross sections is spoiled due to the singularity when the EWIMP mass 
is heavy compared to the weak gauge boson masses. 

1. Properties of EWIMPs 

The EWIMP dark matter is a neutral component of 811(2)^ multiplet(s). In 
this paper, we consider two cases. One is an SU(2)i triplet fermion whose hyper- 
charge is zero. In this case, the EWIMP is accompanied with the SI]{2)l partner, 
a charged Dirac fermion . They are almost degenerate in mass, and the mass 
difference 6m is caused by the electroweak symmetry breaking. If 6m comes from 
the radiative correction of the gauge boson loops pHj . 



where m is the EWIMP mass, a2 is the SU(2)i gauge coupling, miy(m^) is the 
W{Z) gauge boson mass, and cwisw) = cos 9w {sin 9 w), where 9w is the Weinberg 
angle. The gauge interactions of the EWIMP and its SI]{2)l partner are given by 



where e = v47ra and a is the fine structure constant. The mass difference ^mjad 
is induced by the custodial S'U{2)l symmetry breaking in the gauge sector, and 
^mrad — 0.18 GeV if m ^ mw and mz- Effective higher-dimensional operators may 
also generate the mass difference, however it is suppressed by m^/A^, where A is a 
new particle mass scale. 

Another example of the EWIMP dark matter is a neutral component in a pair of 
SU(2)i doublet fermions with the hypercharges ±1/2. After the symmetry breaking, 
two neutral mass eigenstates, x° and x'n: appear. The lightest one is a candidate 
of the EWIMP dark matter. A charged Dirac fermion x~ is also accompanied with 
them, and they are also degenerate in mass in the SU{2)l symmetric limit. The 
mass differences among them are generated by effective operators via the electroweak 
symmetry breaking. Unlike the triplet EWIMP case, the mass difference is 0{m'^/A) 
and it is not strongly suppressed by A. The gauge interactions of the doublet EWIMP 





(1) 
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dark matter and its partners are given by 

+eJ=Ax + TT^'^fxl ■ (3) 

An example of the EWIMP dark matter is the hghtest neutrahno in the MSSM. 
Neutrahnos {i = 1 ■ ■ - 4) are hnear combinations of the supersymmetric partners of 
gauge bosons and Higgs bosons, Bino (B), neutral Wino (W'^) and neutral Higgsinos 
{Hf, i?2)- While those four fields have SU(2)i® U(l)y invariant masses, they are 
mixed with each other via the electroweak symmetry breaking |17j . 



ZaB + Zi^W"^ + Z,sH', + Z^^Hl . (4) 



(5) 



Coefficients Zij are determined by diagonalizing the neutralino mass matrix, 

/ Ml —mz sw cp mz sw Sjs \ 

M2 mz cw cp -mz cw sp 

—mz Sw cp mz cw cp -/i 

\ mz Sw sp -mz cw sp / 

which is written in the {B, W^, H^, H2) basis. Here Mi and M2 are the Bino and 
Wino masses, respectively, and /i is the supersymmetric Higgsino mass. The variable 
tan P is given by the ratio of the vacuum expectation values of two Higgs fields, and 
cp = cos P and sp = sin/5. The lightest neutralino is Wino-like when M2 <^ \fi\,Mi, 
and Higgsino-like when <^ Mi,M2. These two neutralinos have SU(2)j^ charges 
and are candidates of the EWIMP dark matter; the Wino-like neutralino is a triplet 
EWIMP and the Higgsino-like neutralino is a doublet EWIMP. 

Neutralinos are accompanied with charginos x7 = 2), which are linear com- 
binations of charged Wino W~ and charged Higgsino H~ = H^^ + iJ2/? .^'''J- The 
compositions of charginos are determined by diagonalizing the chargino mass matrix, 

\ V2 mw Cp fx J 

which is written in the {W~,H^) basis. 

From the above matrices in Eqs. (jSJ and the mass difference ^mtree between 
the lightest neutralino and chargino at tree level can be calculated. If the LSP is 
Wino-like {mz,M2 <^ Mi, (5mtree is approximately given by 

m'^ 

Smtree ^ J^^w^w si^^ 2/3 , (7) 
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Figure 1: Contour maps of the lightest neutralino mass (solid line) and the mass difference 
between the lightest neutralino and chargino (dashed line) in (M2, /x) planes with tan (3 = A 
(two top figures) and tan /? = 40 (two bottom figures) in the MSSM. M2 = 2Mi is assumed. 
Shaded areas correspond to the Higgsino-like region (l^isp + |-Z^i4p > 0.9). 

which is suppressed by the third power of the high energy scale Mi/x^ as discussed 
before. Since their masses are highly degenerate at tree level, the radiative correction 
to the mass difference in Eq. is also important. 

The mass splitting for the Higgsino-like LSP in a case with mz, |/^| Mi, M2 is 
0(m|/msusY) and given by 

Sm ^ l^/wi^ - sin 2/3) + l^/wi^ + sin 2/3) . (8) 

The second lightest neutralino is also degenerate with the LSP and the chargino in 
mass, because they are in common SU(2)i multiplets. The mass difference 6171^ 
between the LSP and the second lightest neutralino is again (9(m|/msusY), 

Sm, ^ gc^^ + ^4 , (9) 

and this is 2 x 6m when tan/3 ^ 1. 

In Fig.Qand Fig.|2l we show contours of the lightest neutralino mass and the mass 
difference between the neutralino and the lightest chargino in {fi, M2) planes with 
tan /3 = 4, 40. These figures are obtained by diagonalizing the mass matrices in Eqs. 
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Figure 2: Contour maps of the lightest neutralino mass (solid line) and the mass difference 
between the lightest neutralino and chargino (dashed line) in (M2, /x) planes with tan /3 = 4 
(two top figures) and tan/3 = 40 (two bottom figures) in the MSSM. M2 = Mi/ 3 is 
assumed. Lighter shaded areas correspond to the Higgsino-like region (|^i3p + |-^i4p > 
0.9), while the darker shaded areas are the Wino-like one (|Zi2p > 0.9). The radiative 
correction in Eq. is included for depicting the contour of the mass difference in the 
Wino-like regions. 

(0) and (jni) numerically. In Fig. [TJ we assume the GUT relation between the gaugino 
masses, M2 = 2Mi. In this case, the Higgsino-like neutralino (doublet EWIMP) may 
be the dark matter if |yu| ^ Mi. The shaded areas in these figures correspond to the 
Higgsino-like region (|2'i3p + |^i4p > 0.9), and the lightest neutralino is degenerate 
with the lightest chargino, especially at a large mass. 

Fig. 12 is the same plots except that we assume the relation M2 = Mi/3, which 
is predicted in the anomaly mediated supersymmetry breaking scenario [21] . In this 
case the Wino-hke (triplet EWIMP) or the Higgsino-like (doublet EWIMP) dark 
matter may be realized. The lighter shaded bright as the shaded areas in 

Fig.lII) are the Higgsino-like region (|Zi3p+|Zi4p > 0.9), and the darker shaded areas 
are the Wino-like one (|^i2p > 0.9). The Wino-like neutralino is highly degenerate 
with the lightest chargino in mass as expected. 

When the triplet EWIMP mass is around 1.7 TeV, the thermal relic density of 
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Figure 3: Dominant diagram in the Wino- or Higgsino-like neutralino annihilation to 
two photons at one-loop level, when the neutralino is heavy compared to the weak gauge 
bosons. 

the dark matter is consistent with the WMAP data. In the doublet EWIMP case, 
the mass around 1 TeV explains the WMAP data |2S]- However, note that the 
dark matter in the universe may be produced thermally or non-thermally P7] . 
Therefore, we do not assume any scenarios for the dark matter relic density in this 
paper. Instead, we assume the dark matter forms the dark halo in our galaxy with 
the appropriate mass density. 

2. Threshold singularity 

In the calculations of the EWIMP annihilation cross sections, a threshold singu- 
larity appears due to the gauge interactions. For investigating the singularity, let us 
consider the EWIMP annihilation cross section to two photons oi}^^ —>■ 77) as an 
example. The signal line gamma ray flux from the galactic center is proportional to 
~^ 11) ■ This process is a radiative one, and the full one-loop non-relativistic 
cross section in the MSSM context has already been calculated in Ref. [ISj. It is 
found that the cross section is suppressed only by the W boson mass, not by the 
neutralino mass as 



if the neutralino is heavy and almost Wino- or Higgsino-like. The dominant diagram 
is shown in Fig. El 

On the other hand, the cross section must be bounded from above by the unitarity 
limit. 



av < 




(11) 
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Figure 4: Dominant diagram in the Wino- or Higgsino-like neutralino annihilation at 
0(002), which n weak gauge bosons are exchanged. 

Thus, the one-loop cross section exceeds the bound for the extremely heavy neu- 
tralino. It means that the higher-order corrections should be included. The domi- 
nant higher-order contribution comes from the ladder diagrams. The n-th order ) 
ladder diagram, in which n weak gauge bosons are exchanged, is depicted in Fig. |3J 
The corresponding amplitude An of the diagram is roughly given by 



When the neutralino mass m is large enough, the diagrams are enhanced by a factor 
of a2m/my/ for each weak gauge boson exchange. The higher-order loop diagrams 
become more and more important when ^ mw 

Enhancement of ladder diagrams in non-relativistic limits is related to a threshold 
singularity. Recall that a threshold singularity appears in the non-relativistic yU"*'/i~ 
pair annihilation cross section. When the relative velocity v of the muon pair is 
smaller than a, the amplitude of the n-order ladder diagram, in which n photons are 
exchanged between the muon pair, is proportional to ^(a/f)", and the perturbative 
expansion by a breaks down. The internal muons are close to non-relativistic on- 
shell states. The muon and photon propagaters are proportional to and each 
loop integration gives av^ . Thus, the diagrams are enhanced hy a/v for each photon 
exchange. This is because the kinetic energy of muon pair, is smaller than 

the Coulomb potential energy, a^m^, and the wave function of the incident particles 
is deformed from plane waves. We need to systematically resum the ladder diagrams 
or to use the wave function under the Coulomb potential in order to get the precise 
annihilation cross section. 

In the non-relativistic EWIMP pair annihilation, the sub-diagram corresponding 
to the process x^'x^ ~^ X^X~ each ladder diagram is very close to the threshold 




(12) 
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when the mass difference 6m is neghgible. In this case, the spatial momentums for 
EWIMPs and the 811(2)^ partners in the internal lines are regularized by the weak 
gauge boson masses. Their propagaters and the weak gauge boson ones behave as 
m/m^ and l/m"^, respectively, and the loop integration gives a2'm^^/m. Thus, 
the diagrams are enhanced by a2m/mw for each weak gauge boson exchange, as 
shown above. This implies that when a2m^ mw, the weak interaction becomes a 
long-range force and the wave function is significantly modified inside the Yukawa 
potentials induced by the weak gauge boson exchanges in the non-relativistic limit. 
In the following sections, we will introduce a systematic method to derive the annihi- 
lation cross sections in the threshold singularity region by using the non-relativistic 
action. 

The elastic scattering cross section of dark matter with nucleon is important 
for the direct search for dark matter [2H]- If the dark matter is an EWIMP and 
much heavier than the weak gauge boson masses, the one-loop correction to the 
cross section is suppressed only by the weak gauge boson masses and it may be 
dominant over the tree-level contribution [221 • However, the perturbative expansion 
is still reliable, unlike the case of the annihilation cross sections. This is because the 
reduced mass in the EWIMP and nucleon two-body system is not heavy enough for 
non-perturbative corrections to be sizable. 

Ill Two-body State Effective Action 

In this section we derive the effective actions for the non-relativistic two-body 
states including pairs of EWIMPs and the SU(2)l partners. The action is derived 
by following steps, (i) We integrate out all fields except EWIMPs and the SU(2)i 
partners such as the lightest neutralino and chargino in the MSSM action, (ii) The 
non-relativistic action iSnr is obtained by integrating out large momentum modes of 
EWIMPs and the SU(2)l partners. The action includes the effect of the EWIMP 
annihilation as the absorptive parts, (iii) The action 5nr is expanded by the veloci- 
ties of EWIMPs. (iv) At last, we introduce auxiliary fields to the action iSnr, which 
represent two-body states of EWIMPs and the SI]{2)l partners. The two-body state 
effective action S^^^^ is obtained by integrating out all fields except those auxiliary 
fields in the non-relativistic action iSnr. 

In the following, we derive the two-body state effective action for the triplet 
EWIMP (Wino-like neutralino). For the doublet EWIMP (Higgsino-like neutralino). 
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only the final result is presented, because the derivation is essentially the same as 
that of the triplet one. 

1. Integrating out all fields except x° and x~ 

Relevant interactions to the annihilation cross section for the triplet EWIMP 
are the gauge interaction in Eq. (j2)). In the MSSM, the lightest neutralino and 
chargino have other interactions with sfermions and Higgs bosons. However, the 
contributions to the annihilation cross sections are suppressed by the sfermion masses 
or the gaugino-Higgsino mixing. We assume that these contributions are small, and 
neglect them. 

After integrating out the gauge bosons, the effective action for the triplet EWIMP 
and its charged SU(2)l partner becomes 

+ 5int[x°,r], (13) 




where the parameter nic is the mass of the charged SU(2)l partner. The functions 
G^P{x),Gil\x) and CiZJ (x) are the Feynman propagaters of the W , Z bosons and 
photon, respectively. The currents j(^(x) and j^_(x) are defined as 

iwi.^) = X^{x)Y}C{x) , Jx-(^) = X"(a;)7^X~(a;) • (15) 

Here we include the effects of the non-vanishing mass difference between the 
EWIMP and its 811(2)^ partner, 5m{= rric — m), which comes from the electroweak 
symmetry breaking. The non-vanishing Sm gives sizable effects on the annihilation 
cross sections when 6m is not negligible compared with a2mw as will be shown in 
Sec. V. On the other hand, the SU(2)/^(8)U(l)y breaking in the gauge interactions 
gives at most corrections up to 0(m^/A^) to the cross sections, and they can be 
ignored as far as mw <^ A. 

2. Integrating out large momentum modes of x° and x~ 

We now derive the action which describes the non-relativistic motion of the 
EWIMP and its SU(2)i partner. Namely we integrate out the large momentum 



>efr 



d X 



m^)X 
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modes of these particles in «Seff. We divide the fields x° and x into two parts, the 
non-relativistic part and the other, 



(2^ 



<f{p)e 



ipx 



4>%p)e-'^^ (16) 



'[NR] l^'U ^[NR] 

where the (p'^{p) is the Fourier coefficient of the EWIMP field. The region of the 
integration [NR] is defined roughly by [NR] = | P° = ^rn + 0{m\v\'^), p = 

0{mv), \v\ <^ 1}, and [NR] means the complementary set of [NR]. The SI]{2)l 
partner field is also divided into Xt^^ iii the same way. After integrating 

out large momentum modes 5x^ and 6x~ in the action S^s, the non-relativistic 
eff'ective action 5nr is obtained as 



>NR 



d X 



1^ 



X°NR {i^ - m) Xnr + X NR {i^ - mc) Xnr 



+'5int[XNR, Xnr] + <^'5[XnR, Xnr. 



r,0 



[17) 



All effective interactions induced from the integration by 5%° and 6x^ are included 
in SS. Though they are suppressed by the EWIMP mass in comparison with in- 
teractions in iSint, they give leading contributions to the imaginary part (absorptive 
part) of the non-relativistic action. The action is further simplified in the next step. 

3. Non-relativistic expansion of the action iSnr 



Here, we expand the action <Snr by the velocity of the EWIMP. For the expansion, 
it is convenient to use two-components spinor fields (, r] and ^ instead of x^r ^nd 
Xnr- These spinor fields are defined by 



Xnr 



v 



le 



2m 



Xnr 



2m ^ / 



v 



Q—irnt^ _|_ j^^imt 



r,imt( 



le 



-imt 



2m 



(18) 



2m 



The spinor is the charge conjugation of = —ia'^(*, where (T2 is the Pauli 

matrix. Spinors ( and r] annihilate one x° and one x~ , respectively, while ^ creates 



one x^. 
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The non-relativistic action 5nr is systematically expanded by the velocity of the 
two-components spinor fields. The kinetic terms in Eq. ()17|) become 



NR. I kinetic terms 



d X 



Sm H 

2m 



+ e ( + 6m- e 



.(19) 



The interactions in 5int of Eq. (fTTjl are reduced as 



int 



d'^xd^y 



a 



2\x — y\ 

mw\x-y\ 



1 + ^e-^z\^-y\ r^\x)i{y,x') x i\y,x')r^{x) 



026 



{C^(x)C(y,x°) xe^(y,x°)r/(x)+h.c.} 



,(20) 



2\x — y\ 

where we keep terms which dictate the transitions between states with both spin 
and isospin singlet. The EWIMP S-wave state is spin-singlet due to the Majorana 
nature, and only those terms are kept to calculate the EWIMP annihilation cross 
sections at the non-relativistic limit. The first term in the parenthesis describes the 
Coulomb force and the force by one Z-boson exchange between x"*" and ■ The 
second term is for the transition between and x^x~ by one W boson exchange. 

The imaginary part (absorptive part) of the action iSnr is in 5S in Eq. (fT7j) . 
It comes from following box diagrams; (a) transitions from x^x~ to x^xT with 
intermediate W~^W~ , Z'^Z^, 7Z and 77 states, (b) a transition from to x^x^ 
with an intermediate W~^W~ state and (c) a transition from to with an 

intermediate W~^W~ state. These effective interactions are simplified in the non- 
relativistic expansion as 



5S 



d X 



1 



+ + 2s^C^ + sly ] 77"^^ ■ ^^T] 



(21) 



Here, we assume that the EWIMP mass m is much heavier than the weak gauge 
boson masses, m 3> mw,TTiz- The first term of the right-hand side in Eq. ()21|) 
corresponds to the box diagrams (a). Each term in the parenthesis comes from the 
diagrams with intermediate W~^W~ , Z^Z^, 7Z and 77 states, respectively. The 
second and third terms correspond to the diagrams of (b) and (c), respectively. 



4. Two-body state effective action 
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The non-relativistic action for the triplet EWIMP is now given by Eqs. ()19I - I^T|) . 
We now introduce auxihary fields (Tat and ac, which describe the two-body states 
and x^x~ with both spin and isospin singlet, respectively. We thus insert 
identities 



1 



VaNT^Sj^ exp 



VacT>s]j exp 



^ / d^xd^y ac{x,y) \s]j{y,x) -r]\x)^{y,x^) 



, (22) 



and their conjugates into the partition function described by the non-relativistic 
action. After integrating out t], ^, (, sn, sc and their conjugates, the two-body state 
effective action S^^^^ is obtained as 

= J d'^xd^r f) I (^td^o + ^ + - V(fO + 2ir(5(f)| <l>(x, r) , (23) 

where the argument x denotes the center of mass coordinate in the two-body system 
and r is the relative coordinate. The two- components two-body state field $(x, r) is 
defined by 



$(x, r) 



f)c[x,r) 
')N{x,r) 



-1 / (^c[x,r) 
aNix,f) 



(24) 



The components, (j)^ and (j)c, describe pairs of the EWIMPs and the 811(2)^ partners, 
respectively. These can be also written as 



1 



(j)c{x,r) = —=^^(x — r/2,x^)r]{x + r/2,x^) 
v2 



(25) 



^N{x,f) = ^C\x-f/2,x°)C{x + f/2,x'^) 



(26) 



by using (, t] and ^ fields. The difference between the normalizations of (pc in 
Eq. and 0Ar in Eq. ()26|) comes from the fact that (p^ describes a pair of the 
identical Majorana fermions. 

The electroweak potential V(r) in Eq. ()23|1 is given by 



Vfr) 



a 



-mzT 



25m Ci2Cw 



-V2a2- 



\ 



-V2a2- 



(27) 
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The absorptive (imaginary) part T is decomposed as T 
where each component is 



w+w- 



■ 77 



w+w- 



■yZO 




Z"ZO 



■ 77 




(28) 



The two-body state effective action Eq. 



is the final result of this section. 



In the case of the doublet EWIMP, the derivation of the effective action is parallel 
to that of the triplet one. The EWIMP is accompanied with the neutral S\J{2)l 
partner in addition to the charged one. The two-body effective action becomes 3x3 
matrix form, and the two-body state field $(x, r) has three components, 



$(x, r) 



4>c{x,r) 



\ 



0Ar(x,f) 



(29) 



where (pN, (pc and (pjy^ describe pairs of the EWIMPs, the charged and the neutral 
partners, respectively. The electroweak potential is 



/ a 02(1 -2c^y 

2dm 

r 4c 



V(r) 



2 \2 ^-mzT 



^2, 



-mw'T 



^2, 



V2i 



026 



-mwi' 



72, 



4r 

0^26 



4r 




4r 



Q!2e 



-mzT 



Ac^r 



-razT 



\ 



4r 



4c^r 



25m 



N 



(30) 



where brnfq is the mass difference between the EWIMP and its neutral SU(2)/^ part- 
ner. The absorptive part F in this case is given by F = F^y+vy- +F^020 + F^^o + F^^, 
where 

/ 8 72 



w+w~ 



nan 



nan 



V2 4 
\^v^ 4 
/ 4(1 



ZOZO 



■yZO 



naa^ 
2c^m2 



4 
4 



72(l-2s^)2 



(il- 



\ 



v^(l - 2.2^)2 
2 
2 



v^(l 



2^^)A 



2 
2 



'isyf o\ 

oy 



■ 77 



na 



J 

A o^ 

000 

oy 



.(31) 
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IV Optical Theorem and Annihilation Cross Sections 



We now derive the EWIMP pair annihilation cross sections by using the two- 
body state effective actions in Eq. (j^^ . Using the formula, we also show that the 
one-loop cross section of x^x^ 27 is reproduced in the perturbative expression 
while the cross section in a limit of m — > 00 satisfies the unitarity bound. 

1. Annihilation cross section formula 

Due to the optical theorem, the total EWIMP pair annihilation cross section is 
written by the imaginary part of the amplitude as 



Vs^ - ^m^s a = $5 [Mii] , (32) 

where y/s is the center of mass energy, and J^a is the invariant amplitude of the 
process, ~^ X^X^- When the incident EWIMPs are highly non-relativistic, the 
relevant initial state for the pair annihilation is only the spin-singlet S-wave state 
as mentioned before. We thus project the above equation (jH^ to the ^Sq state by 
using the projection operator Jd^Pdk\P,k, ^So;N){P,k, -'^5'o;A|, 

|P, k, ^So; N) = / dilk \ aiiP/2 + k) al(P/2 - k) 

4V7r J L 

-al(P/2 + k) a\{P/2 - k)] |0) , (33) 

where P is the total momentum and k is the relative momentum in the EWIMP two- 
body system. The variable k is k = \k\. The state vector \P,k,^ So) is normalized 
as {P',k', ^So;N\P,k, ^So;N) = 6{P - P')6{k - k'). The operator a}^{k) creates 
one x° with the momentum k and the spin ±1/2, and satisfies the anti-commutation 
relation {a^(p), a^,(p')} = 6{p-p')6ss'- 

If the x~^x~ annihilation cross section is considered in a spin singlet and S'-wave 
system, the projection operator J d^Pdk\P, k, ^So;C){P,k,^So;C\ is used instead 
of one in Eq. (jHHj) . where 

\P,k,^So;C) = —^[dnk\ bl{P/2 + k) dl{P/2-k) 
2v27r J L 

-bl{P/2 + k) dl{P/2 - k)] |0) . (34) 

The operators &±(A;) and dl_(k) create one x~ ^ind one with the momentum 
k and the spin ±1/2, respectively. They satisfy the anti-commutation relations 
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After the projection, the S'-wave cross sections {i = 1,2) become 



is) _ 327r5 ^ 

— Cj ^ 



(35) 



( s) 

where v is the relative velocity between incident particles. The cross section is 
for annihilation {i = 2) and X'^X~ annihilation {i = 1). The coefficient Cj is 

given by ci = 1 and C2 = 2. Ai[^\v) is the invariant amplitude from the ^5*0 state 
to the ^5*0 state. We use the non-relativistic approximation as s ~ 4m^ + m^f^ for 
deriving the cross sections. 

The invariant amplitude Ai\^\v) is obtained from the S-wave Green function of 
the two-body state field From the two-body state effective action in Eq. (j23p . the 
Schwinger-Dyson equation (equation of motion for the Green function) is derived as 

(zd^o + + Z? _ v(r) + tV^] {0\T^x, r)$t(y, r')\0) 

= i5^^\x-y)5^'^\f-f^) . (36) 

The potential V(r) and the absorptive part F are defined in Eqs. ipTj) and (j^H|l for 
the triplet case, and in Eqs. (pUjl and (PT|) for the doublet case. Since the potential 
depends on only r, the Green function can be expanded by the Legendre Polynomials 
as 

p or _|_ 1 

__e-^(-i/)^_^P^(cos7)(-0G(^'')(r,r') , (37) 

where 7 is the angle between r and r*. The variable E in the superscript is the 
internal energy of the two-body state {E = — P^/Am). The equation for the 
function G^^'''\r,r') is then given by 

.E-±f^r-'Sl±^^ V(.) - ,rM) G<-.(. ^ ^ . (38) 

The invariant amplitude }A\ \v) is written by the (i, i)-component of the S'-wave 
Green function G-f as 

Mf\v) = ^, lim [e / r^c/r r''rfr'jo(A;r)jo(A;r')Gjf'°^(r,r') , (39) 

where k = mv/2 and jo{x) is the zero-th order of the spherical Bessel function. By 
using this expression, we obtain the formula for the total S-wave annihilation cross 
section al^\ 

(Tp^f = Cj— lim [E ) / rdr r'dr' sm(kr) sinikr')'^ G^f'^Vr, r') .(40) 

k^ E^k^/m \ rn J Jq L J 
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In general, the function G ' (r, r') can not be solved analytically, therefore we need 
to solve the Schwinger-Dyson equation (PH|) numerically to obtain the cross sections. 

2. Solving the Schroedinger equation 

When a function, g(r, r') = rr'G*^^'°''(r, r'), is defined, Eq. ()38p becomes the 
standard Schroedinger equation in one dimension, 

- + f^(^) - g(r, r) - Eg(r, r') = 5(r - r') . (41) 

In the following we expand the solution g(r, r') perturbatively by the absorptive part 

r. 

At the leading order, the solution go(r, r') satisfies the equation, 
1 

— go(r, r') + V(r)go(r, r') - ^go(r, r') = 5(r - r') . (42) 

The boundary conditions of the equation are determined by the following two re- 
quirements; (i) The Green function G'''^'°^(r, r') must be finite for any r and r' . (ii) 
The Green function has only an out-going wave at |r — r'| ^ oo. Then the solution 
go (r,r') is obtained as 

go(r, r') = mg>(r)g^(r')6'(r — r') + mg<(r)g^(r')6'(r' — r) , (43) 

where g>(<)(?") is the solution of the homogeneous equation, 
1 

- -^g>«) W + vWsxo W = ^g>«) W ■ (44) 

Here, g>(<) is given by a 2x2 (3x3) matrix in the triplet (doublet) EWIMP case. The 
solution g<(r) satisfies the boundary conditions, (i) g<(0) = and (ii) g<(0) = 1, 
while the conditions for g>(r) are (i) g>(0) = 1 and (ii) g>(r) has only an out-going 
wave at r ^ oo. 

At the first order of F, the solution of Eq. ((H)), gi(T', r'), is simply given by the 
leading-order solution go(T', t') as 

gi(r,r') = rfr"go(r,r") ("^r^^ go(r",r') = ^g>(r) V g>(r') . (45) 

The S'-wave cross sections in Eq. fl4()j) are proportional to the imaginary part of 
the Green function G^-f'''''. The imaginary part is related to not only the annihila- 
tion cross section but also the elastic cross section of the process, oi' 
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X^X ~^ X^X ■ After extracting the contribution of the annihilation processes from 



53 



G 



{E,0) 



, the total S'-wave annihilation cross section is obtained as 



a] V = c 



fc2 E^k^/r 
oo 

Aia = I rfrsin(/cr) [g>(r)]^^ . (46) 







It is found that only the asymptotic behavior of the function g>(r) is relevant 
in the calculation of the cross sections. The function g>(r) has only an out-going 
wave at r ^ oo as stated above. When the SU(2)l partners do not appear in the 
asymptotic state {E < 26m), the function g>(r) should behave as 

[g>W],, = 5,2 d,,iE) e^^^ . (47) 



Thus, the total EWIMP pair annihilation cross section is given in a simple form, 

2 J2 d2aimvyA) dl,{mv^/A) . (48) 



(To V 



If the potential term in Eq. ()42|1 is neglected (in other word, the long-distance 
effects are negligible), the coefficient d2a{E) is given by 62a and the cross section is 

(S) 

simply given by cig v = 2T22, which is consistent with the tree-level cross section as 
expected. Below we omit the subscript 2 for the EWIMP annihilation cross section 
for simplicity. 

3. One-loop result in x^x*^ 77 process 

For demonstration of the validity of Eq. P8|l , we show that the EWIMP annihi- 
lation cross section to two photons agrees with the one-loop cross section when the 
potential term V is treated perturbatively. 

From Eq. ()48|). the annihilation cross section to two photons is given by 

a^'^^^ = 2[rg,^ \d2iimvyA)\\ (49) 

where the partial absorptive part to two photons F^^ is given in Eq. ()28|) for the 
triplet EWIMP case and in Eq. (|3ip for the doublet EWIMP case, respectively. Here 
we show only the result in the triplet EWIMP case for simplicity. The coefficient ^21 
is obtained by solving the Schroedinger equation ()44|) . When we expand the solution 
by the potential term V(r), [g>(r)]2i is obtained at leading order as 

[g>(r)]2, = d2,iE) ^ -mv^2a2 (^q) 

mw + V2mdm 
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Here, we take E < 26m and m ^ niw- Thus, the annihilation cross section to two 
photons is given as 

-'^'H. -^^fl + ^/^V'. (51) 



This agrees with the result obtained in the full one-loop calculation in a heavy 
EWIMP mass limit |19|, including the correction due to the non- vanishing 6m. 

4. Cross section for x^x^ ^ 77 in a limit of m 00 

When the EWIMP mass m is much heavier than the weak gauge boson masses, 
we can not deal with the potential term perturbatively. The equation (jl^ can be 
solved analytically in the limit of m 00, and the qualitative behavior of the cross 
sections can be discussed. This is because the weak gauge boson masses mw and 
mz and the mass difference 6m can be neglected in this limit. The Schroedinger 
equation for g> (r) becomes 

- --f^sAr) + -Ug>(r) = Eg>(r) , (52) 
m dr^ r 

where U is the coefficient matrix for the electroweak potential, and defined by U = 
[rV(r)],.^o- Because all forces in the potential become Coulomb-type, the solution 
of above equation is determined by using the confiuent geometric function as 

V 

^•('■' = ^ ^yi) ^-^^A (-^'^O • (53) 

The function r(z) is the Euler's Gamma function and W^^^ii^z) is the Whittaker 
function. The matrix O is the diagonalization matrix, 0^\J0 = diag.(Ai, ■ ■ ■ )' 
Xi is the eigenvalue of the matrix U. 

The EWIMP annihilation cross section to two photons is then derived as 

2 



2.8 X IQ-^vm^ (Triplet EWIMP) 
3.2 X IQ-'^/vm^ (Doublet EWIMP) ' 
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in the limit of m — oo. Here, we use the approximation f ^ 1 to derive the last 
equation. The cross section behaves as av ~ l/(fm^), and satisfies the unitarity 
condition in Eq. (fTTj) as expected. 

V Annihilation Cross Sections and Zero-Energy Resonances 

When the EWIMP mass is heavy enough so that the the effects of the long- 
distance force by the electroweak potential cannot be ignored but not heavy enough 
to take the limit of m — oo, we have to solve the Schroedinger equation ()44|) nu- 
merically for the precise annihilation cross sections. In this section we show some 
numerical results, and discuss the behaviors using a toy model. 

1. Numerical result and zero-energy resonances 

First, we show numerical results for the triplet and doublet EWIMP annihilation 
cross sections to 77 and W^W~ . In Fig. El the cross sections with some fixed mass 
differences {5m = 0.1,1 GeV) are shown as functions of m. We set to the mass 
difference 5mjv to be 6171^ = 25m for the doublet EWIMP. In this calculation, the 
relative velocity of the incident EWIMPs v/c is taken to be 10"'^, which is the typical 
velocity of dark matter in the galactic halo. We also show the leading-order cross 
sections in perturbation as dotted lines. When the EWIMP mass m is around 100 
GeV, the cross sections to 77 and W^W~ are almost the same as the perturbative 
ones. However, when m is large enough (m^ 0.5 TeV for the triplet EWIMP and 
1.5 TeV for the doublet EWIMP), the cross sections are significantly enhanced 
and have resonance behaviors. 

The qualitative behavior of the cross sections, especially around the resonances, 
may be understood by using a toy model, in which the electroweak potential is 
approximated by a well potential. Here we discuss the triplet EWIMP case, and the 
extension to the doublet case is straightforward. Taking cw = 1 and 5m = for 
simplicity, the electroweak potential is approximated as 



where bi and 62 are numerical constants. By comparing the annihilation cross sec- 
tions to 77 in this potential and in the perturbative calculation Eq. 1)511) for small 
m, we find bi = 8/9 and 62 = 2/3. Under this potential, the two-body states of 




(55) 
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Figure 5: Annihilation cross sections (av) to 77 and l^+VF" when dm = 0.1, 1 GeV for 
both the triplet and the doublet EWIMPs. Here, v/c = 10~^. The leading-order cross 
sections in perturbation are also shown for Sm = (broken lines). 

2x° and x'^x~ have attractive and repulsive states, whose potential energies are 
A± = (Vii ± VVfi + 4V?2)/2 with {ij = 1,2) elements in V. The attractive 
state is cos^ 0jv — sin^ (f)c with tan^ 9 — — A_/A+. 

By virtue of the approximation, the pair annihilation cross sections for the triplet 
EWIMP are obtained analj^ically. 



(av) 



w+w- 



7ra| 
9m^ 



21(^22) + 21^22!^^ , (av) 



^21 = I V2 
d'2.2 = I 

+2 



cos 



v2 + 2vl 



sin [vc^Jv^jvl + 2^ 



-1 



cos PcV ^ / ^c ^ 1 



cos [Pc\/v'^/vl + 2 



23^ sin (peVt'V^'c - 1) 
Vv7^;2 + 2) 



t'2 + 2?;2 



sm Pc 



-1 



cos 



Vc\Jv^lvl - 1^ 



-1 



(56) 



where Pc ^-nd I'c are defined by Pc = ■\/2a2m/mw and = •\/32Q;2mvK/9m, respec- 
tively. 

If the kinetic energy of the EWIMP pair is much larger than the potential energy 
{v » Vc) or the electroweak potential is point-hke (pc ^ 1), the cross sections 
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coincide with the results in leading-order calculation in perturbation as expected. 
However, when v is much smaller than Vc, such as in the case of dark matter in the 
current universe, the d2i and ^22 become 



d2i — cos a/2pc —V2 [cosh Pc] 

^^22 



cos\/2pc + 2 [coshpc] . (57) 

Therefore, when V2pc — {2n — l)7r/2 (n = 1, 2, ■ ■ ■ ), the coefficients ^21 and ^22 are 
enhanced by several orders of magnitude. As the result, the cross sections show the 
resonance features as in Fig. These resonances are called zero-energy resonances 
pUj . because the condition \/2pc — {2n — l)vr/2 is nothing but existence of a bound 
state with zero binding energy. The bound states consist of mixtures of pairs of the 
EWIMPs and the charged partners. 

In Fig. the first resonance in small 5m (that is Sm ~ 0.1 GeV) appears at 
m ~ 2 TeV. On the other hand, the well potential model predicts the first resonance 
at m ~ 1.8 TeV. Thus, the model describes the behavior around the first resonance 
well. Also, notice that the cross section to two photons, which is induced by the 
one-loop diagrams in perturbation, is suppressed only by j(x^ compared with that 
to W~^W~ for ~ 1- This behavior is also seen in Fig. El 

On the resonance, the coefficients ^21 and ^22 in Eq- behave as 

d21 ^ t— , (322 - ^ , (58) 

V V 

for the small relative velocity v <^ Vc- Thus the cross sections av are proportional 
to v""^. However, this is not a signature for breakdown of the unitarity condition in 
Eq. (fTT|. We find from study in the one-fiavor system under the well potential V 
that av would be saturated by the finite width for the two-body system (that is, the 
absorptive part) F when v <^ mVT, and the unitarity is not broken. 

So far we have ignored the mass difference 6m in Eq. ()55p . When 6m is not negli- 
gible compared with ~ a2mw, the potential energy for the attractive state is reduced 
as |A_| ~ lQa2mw/9 — A6m/3 for 6m^ a2miv and |A_| ~ 128{a2mw)'^ / {816m) for 
6m^ a2mw, and the attractive force becomes weaker. Also, the component of the 
EWIMP pair in the attractive state becomes smaller. Then, the zero-energy reso- 
nances move to heavier m for larger 6m, and the EWIMP annihilation cross sections 
around the resonances become smaller, as in Fig. El This effect is more significant in 
the doublet EWIMP case than in the triplet one, since the potential energy by the 
weak gauge boson exchanges is smaller. 
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Now we have discussed importance of the zero-energy resonances for the EWIMP 
annihilation cross sections. One might consider that the thermal relic EWIMP abun- 
dance is also modified due to the zero-energy resonances. The answer to the question 
is "No". The typical relative velocity is given by v ~ 1/3 at the freeze out temper- 
ature T ~ m/20. Note that the critical velocity Vc for the triplet EWIMP is given 
by Vf. ~ 2^2 ~ 0.07 at the first zero-energy resonance where \/2pc = 7r/2. The 
critical velocity for the doublet EWIMP is smaller than that of the triplet one. Thus 
the relative velocity at the freeze-out temperature is not small enough to affect the 
thermal relic abundance. 

One might have another question relevant to the annihilation in the universe 
at very low temperature. Since the cross sections arc proportional to f on the 
zero-energy resonance, the EWIMP annihilation rate (the annihilation cross section 
multiplied by the EWIMP number density Uew) becomes much larger after the 
freezing-out phenomenon stated above. If the annihilation rate becomes larger than 
the Hubble constant H at low temperature, the EWIMP begins to annihilate again 
(re-coupling). The usual calculation of dark matter abundance might be changed. 

For studying the possibility, we use the previous toy model for the triplet EWIMP. 
The thermal-averaged annihilation cross section {av) on the first zero-energy reso- 
nance at temperature T is derived as 



where g ~ 7 x 10 ^. It is found that the condition of the 're-coupling', {av)nEw > H, 
attributes to the inequality. 



where pEw(= mnEw) is the EWIMP mass density and s is the entropy of the uni- 
verse. The lower bound is larger than the cosmological observation for dark matter 
{p/s ~ 7 X 10~^ GeV). Thus the re-coupling does not occur in our universe. 

2. Fitting functions 

The numerical calculation of the EWIMP annihilation cross sections takes huge 
CPU power. It is thus convenient to derive the fitting functions which reproduce the 
numerical results with enough precision. We performed two parameter fitting of the 
annihilation cross sections for the doublet and triplet EWIMPs, the EWIMP mass 




(in the unit of GeV"^) , 



(59) 



Pew 



> 7 X 10"^ [GeV] , 



(60) 



s 
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m and the mass difference between tfie EWIMP and its cfiarged partner 5m. Tfie 
velocity dependence of the cross sections is weak except for very narrow regions in the 
vicinity of the resonances and ignored. For the doublet EWIMP, the cross sections 
further depend on the mass difference between the EWIMP and its neutral partner 
SniN- We fix Suin = 25m in the derivation of the fitting functions for simplicity. 
This relation is valid in the MSSM when tan j3 is large. 
We use the fitting function, 

av{m,5m)\^^ = 10"^*[cm^ s"^] x 

in the range 0.2 TeV < m < 10 TeV and 0.1 GeV < 5m < 1 GeV. By performing the 
least-square method between the fitting function and the numerical integration, the 
coefficients a^- in each process are obtained. The numerical values of the coefficients 
aij for the annihilation cross sections to 77 and W^W~ in the triplet EWIMP case 
and 77, W~^W~ and Z^Z^ in the doublet EWIMP case are given in Table 1. Other 
annihilation cross sections are evaluated from the cross sections to two photons as 

in the triplet EWIMP case, and 

avix'x'^^Z') = 

in the doublet EWIMP case. This is because these processes are induced via the 
transition of the EWIMP pair to the charged partner pair. 

3. Annihilation cross sections in the MSSM parameters 

By using above ffiting functions, we scan the annihilation cross sections on the 
MSSM parameters. The input parameters are the Bino mass Mi, the Wino mass M2, 
the Higgsino mass /j,, and tan/?. We assume that the other super symmetric scalar 
particles, sleptons and squarks, are heavy enough, and neglect the contributions to 
the cross sections. Furthermore we consider two situations, in which relations be- 
tween Ml and M2 are different. The first one is the GUT relation Mi = O.5M2, 
which is frequently assumed in the minimal supergravity scenario. Another one is 



U,.7=0 



m \ ' 



1 TeV. 



5m V 
1 Gev) 



(61) 



= ^^(xV - 77) X , 

a 

= <T.(xV - 77) X ^ , (62) 
'7v{x"x 77 X \ ^ '^^ , 63 
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(Triplet EWIMP, to 77) 





j = 


j = 1 


j = 2 


j = 3 


j = 4 


j = 5 


j = 6 


i = 


2.56521E-1 


2.88649E-2 


1.24874E-1 


-4.97529E-1 


8.48760E-1 


-6.97638E-1 


2.23588E-1 


i = 1 


-1.10191E-1 


8.27605E-1 


-3.13417E+0 


7.52692E+0 


-1.02429E+1 


7.27770E+0 


-2.09385E+0 


i = 2 


-1.59340E-2 


-4.83605E-1 


2.29288E+0 


-5.90629E+0 


8.31712E+0 


-6.03157E+0 


1.75985E+0 


i = 3 


-8.78312E-3 


2.29927E-1 


-1.11637E+0 


2.82055E+0 


-3.86927E+0 


2.73075E+0 


-7.76046E-1 


i = 4 


4.66594E-3 


-5.84155E-2 


2.72064E-1 


-6.65170E-1 


8.86046E-1 


-6.08312E-1 


1.68406E-1 


i = 5 


-5.44363E-4 


6.00290E-3 


-2.73828E-2 


6.56610E-2 


-8.58639E-2 


5.79081E-2 


-1.57571E-2 


i = 6 


1.95407E-5 


-2.10046E-4 


9.48764E-4 


-2.24620E-3 


2.89895E-3 


-1.92989E-3 


5.18540E-4 


(Triplet EWIMP, to W+W) 




j = 


j = 1 


j =2 


j =3 


j = 4 


j = 5 


j = 6 


i = 


-3.36985E-3 


1.23215E-2 


-4.08164E-2 


5.96407E-2 


-2.63576E-2 


-1.42661E-2 


1.13721E-2 


i = 1 


5.61768E-2 


-6.02547E-2 


1.31888E-1 


-3.54816E-2 


-2.97169E-1 


4.06116E-1 


-1.57962E-1 


i = 2 


-4.18499E-2 


8.09353E-2 


-1.05570E-1 


-1.60918E-1 


6.66742E-1 


-7.10708E-1 


2.52266E-1 


i = 3 


6.54865E-3 


-4.19556E-4 


-8.61235E-2 


3.59611E-1 


-6.42137E-1 


5.31727E-1 


-1.66897E-1 


i = 4 


1.43330E-4 


-6.78893E-3 


4.15698E-2 


-1.20926E-1 


1.83077E-1 


-1.37981E-1 


4.08085E-2 


i = 5 


-8.33119E-5 


1.03757E-3 


-5.22051E-3 


1.36654E-2 


-1.93028E-2 


1.38678E-2 


-3.96380E-3 


i = 6 


4.08625E-6 


-4.38210E-5 


2.04733E-4 


-5.07069E-4 


6.86353E-4 


-4.77300E-4 


1.33068E-4 


(Doublet EWIMP, to 77) 


aij 


j = 


j = 1 


j = 2 


j = 3 


j = 4 


j = 5 


j = 6 


i = 


1.04998E+0 


-1.92220E-1 


2.60536E+0 


-8.25766E+0 


1.24256E+1 


-9.11921E+0 


2.61707E+0 


i = 1 


-4.17041E-1 


4.24814E+0 


-1.96145E+1 


5.23145E+1 


-7.39832E+1 


5.28533E+1 


-1.49773E+1 


i = 2 


1.03008E-1 


-1.93097E+0 


1.05626E+1 


-2.97310E+1 


4.29137E+1 


-3.10093E+1 


8.82103E+0 


i = 3 


-2.42293E-2 


4.15616E-1 


-2.28085E+0 


6.20020E+0 


-8.49779E+0 


5.80186E+0 


-1.52947E+0 


i = 4 


2.24803E-3 


-3.43865E-2 


1.66100E-1 


-3.65744E-1 


3.37953E-1 


-9.30267E-2 


-2.82121E-2 


i = 5 


-6.51687E-5 


6.85671E-4 


-2.51709E-4 


-1.22816E-2 


4.22829E-2 


-5.14471E-2 


2.24206E-2 


i = 6 


-1.47604E-7 


1.74390E-5 


-2.40349E-4 


1.19010E-3 


-2.71981E-3 


2.87046E-3 


-1.16460E-3 


(Doublet EWIMP, to W+W) 


aij 


j = 


j = 1 


j = 2 


j = 3 


j = 4 


j = 5 


j = 6 


i = 


-6.51180E-3 


2.79481E-2 


-5.26301E-2 


4.49224E-2 


9.13789E-3 


-4.38570E-2 


2.11067E-2 


i = 1 


1.71330E-1 


-1.78724E-1 


5.13363E-1 


-8.98539E-1 


8.56009E-1 


-3.75230E-1 


4.87411E-2 


i = 2 


-5.42160E-2 


1.85541E-1 


-5.40450E-1 


9.57081E-1 


-8.89508E-1 


3.57243E-1 


-3.19227E-2 


i = 3 


2.44545E-3 


-9.05849E-4 


-3.64007E-2 


1.57625E-1 


-3.13715E-1 


2.94701E-1 


-1.03643E-1 


i = 4 


2.49469E-4 


-3.59210E-3 


2.13727E-2 


-6.52043E-2 


1.08466E-1 


-9.07532E-2 


2.95480E-2 


i = 5 


-2.49263E-5 


3.33127E-4 


-2.02985E-3 


6.39148E-3 


-1.07614E-2 


9.00477E-3 


-2.92073E-3 


i = 6 


6.51335E-7 


-9.63636E-6 


6.24283E-5 


-2.04190E-4 


3.50050E-4 


-2.95151E-4 


9.60757E-5 


(Doublet EWIMP, to Z'^Z°) 


aij 


j = 


j = 1 


j = 2 


j = 3 


j = 4 


j = 5 


j = 6 


i = 


-4.38874E-5 


-9.72321E-3 


7.53513E-2 


-2.21045E-1 


3.21256E-1 


-2.32073E-1 


6.64397E-2 


i = 1 


1.44964E-1 


8.69220E-2 


-5.06389E-1 


1.34705E+0 


-1.86944E+0 


1.31757E+0 


-3.72003E-1 


i = 2 


-6.16430E-3 


-1.64703E-1 


8.62852E-1 


-2.17291E+0 


2.92730E+0 


-2.02460E+0 


5.63888E-1 


i = 3 


-1.36278E-2 


1.29107E-1 


-5.70725E-1 


1.35708E+0 


-1.77804E+0 


1.20897E+0 


-3.32609E-1 


i = 4 


2.05191E-3 


-1.94524E-2 


8.76460E-2 


-2.13319E-1 


2.86459E-1 


-1.99525E-1 


5.61327E-2 


i = 5 


-1.06900E-4 


1.07990E-3 


-5.08739E-3 


1.28877E-2 


-1.79673E-2 


1.29591E-2 


-3.76293E-3 


i = 6 


1.93736E-6 


-2.09370E-5 


1.03637E-4 


-2.74183E-4 


3.97811E-4 


-2.97743E-4 


8.93856E-5 



Table 1: Coefficients of tlie fitting function Eq. for the EWIMP annihilation cross 
sections to 77 and W~^W~ in the triplet case, and to 77, W~^W~ and Z^Z^ in the doublet 
case. 
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Ml = 3M2, which is predicted in the anomaly mediated supersymmetry breaking 
scenario. The hghtest neutrahno mass m and the mass difference between the neu- 
trahno and the chargino 6m in these two cases are given in Fig. ^ and Fig. |21in Sec. 
II. 

The non-perturbative effects due to the existence of the resonances are important 
for the calculation of the annihilation cross sections when the neutralino is Wino- or 
Higgsino-like and heavy enough. On the other hand, the leading-order calculation 
is precise enough in other regions of the parameters. We thus match the fitting 
functions derived in the previous section to the leading-order cross sections in per- 
turbation at m = 250 GeV and 6m = 1 GeV with |^i2p > 0.9 for the Wino-like and 
|Zi3p + |2'i4p > 0.9 for the Higgsino-like neutralino. Here, Zij is defined in Eq. (^j). 

The results are shown in Fig. |H1 and Fig. [71 which are contour plots of the cross 
sections to W^W~ and 77. In Fig. [HI the relation Mi = O.5M2 is imposed. The 
regions |yu| ^ O.5M2 correspond to the Higgsino-like regions. As in Fig. ^ the mass 
difference between the neutralino and the chargino is larger than 1 GeV in most 
of the parameter space. As the result, the first resonance appears at ~10 TeV, 
which is out of the range of these figures. In the regions with O.5M2 ^ the LSP 
is the Bino-like neutralino. 

The relation Mi = 3M2 is imposed in Fig. [7| In addition to the Higgsino-like 
regions, which appears in the same place as ones in Fig. IHl the Wino-like regions 
also appear in the regions M2 ^ The first zero-energy resonance of the Wino-like 
neutralino is shown up at M2 ~ 2 TeV. 

We also show the ratios of the annihilation cross sections including the effects 
of the resonances and the leading-order ones in perturbation in Fig. |H1 and Fig. El 
All parameters for depicting these figures are the same as the those for Fig. El and 
Fig. [71 Huge enhancements are found in the vicinities of the zero-energy resonances. 
Furthermore, the cross sections are enhanced by factors even for m^500GeV when 
the lightest neutralino is Wino-like. Thus, the non-perturbative effect is important 
to obtain the precise annihilation cross sections. The ratios in the Higgsino-like 
regions in Fig. El are larger than that in Fig. |S1 This is because the mass difference 
6m is smaller in Fig. El than in Fig. [HI The non-perturbative effects should also be 
included for the Higgsino-like neutralino with m ^ 500 GeV when the accurate cross 
sections are required. 
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Figure 6: Contour maps of the neutralino annihilation cross sections (av) to W+PF" 
(four top figures) and 77 (four bottom figures) in (M2,//) planes with tan/3 = 4,40 in the 
MSSM. Ml = O.5M2 is assumed. 
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Figure 7: Contour maps of the neutralino annihilation cross sections {ov) to W~^W 
(four top figures) and 77 (four bottom figures) in (M2,//) planes with tan/3 = 4,40 in the 
MSSM. Ml = 3M2 is assumed. 
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Figure 8: Ratios of the cross sections including the non-perturbative effects and the 
leading-order ones in perturbation, (r/fjieading) in the MSSM. Figures are shown as contour 
maps in (M2,/x) planes. All parameters needed for depicting the lines are the same as 
those in Fig. El 
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Figure 9: Ratios of the cross sections including the non-perturbative effects and the 
leading-order ones in perturbation, (r/fjieading) in the MSSM. Figures are shown as contour 
maps in (M2,/x) planes. All parameters needed for depicting the lines are the same as 
those in Fig. [Tj 
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VI Gamma Ray Flux from the Galactic Center 

The enhancement of the dark matter annihilation cross sections has significant 
implication for indirect searches for dark matter using cosmic rays. In this section, 
we discuss the searches for gamma rays resulting from the dark matter annihilation 
in the galactic center and the future prospects. 

1. Flux formula 

The spectrum of gamma rays from the dark matter annihilation consists of two 
components. One is the line gamma rays and the other is the continuum gamma 
rays. The line gamma rays are produced by the radiative processes such as the dark 
matter annihilations to 77 and Z7 |19j. Since the dark matter is non-relativistic 
in the galactic halo, the resulting spectrum is monochromatic. The signal is robust 
for the dark matter search, because the diffused gamma-ray background from the 
astrophysical sources has a continuous energy spectrum. 

The continuum gamma ray signal come from jets in the dark matter annihilation. 
For example, the dark matter annihilates to W bosons, the W bosons fragment into 
TT mesons and finally vr" mesons decay to 77. The energy spectrum thus becomes 
continuous. The observation of the continuum gamma rays may also constrain the 
properties of dark matter if the astrophysical background is understood well. 

The gamma ray fiux from the dark matter annihilation, J-'-y{E), is given by 



where 6 is the angle between the direction of the galactic center and that of ob- 
servation. The function Nj^\e) is the number of photons with energy E in the 
fragmentation of the final state /, and (av) is the dark matter annihilation cross 
section averaged with respect to the velocity distribution function. The density p in 
the integrand is the dark matter mass density profile in our galaxy. 

After integrating out Eq. by the solid angle with the appropriate angular 
resolution of the detector, we obtain 
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Figure 10: Fragmentation functions of W^W~ (the left figure) and Z^Z^ (the right 
figure). The cross points are the simulation data by HERWIG. The fitting functions are 
also depicted in these figures as solid lines. 

The angular resolution of the detector, AQ, is taken to be 10~^ in this paper, which 
is a typical value for the ACT detectors. The information of the dark matter density 
profile is included in a dimensionless function. 



We need three quantities for evaluating the gamma ray fluxes; the dark matter 
annihilation cross sections to final states / ( (crv) f) , the fragmentation functions 
(dNj/dE), and the mass density profile of the dark matter (p). The cross sections 
are obtained in the previous section. We discuss the fragmentation functions and 
the mass density profile in the following subsections. 

2. Fragmentation function 

Since we focus on the SU(2)i non-singlet dark matter, such as the Wino- or 
Higgsino-like neutralino in the MSSM, the continuum gamma ray signal come mainly 
from the the dark matter annihilation modes into W^W~ and Z^Z^. We thus need 
two fragmentation functions, N^+^^ and N^ozo- We simulate the photon spectrums 
from these weak gauge boson states by the HERWIG Monte Carlo code [3T]. We 
derive fitting functions from the simulated fragmentation functions for m = 200 GeV 
and m = 5 TeV by introducing the scaling variable x = E/m. 

In Fig. ^1 the HERWIG simulation data for the fragmentation functions are 
shown. In these figures, the fitting functions are also depicted as solid lines. The 




(66) 
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Figure 11: Dependence of J in Eq. H66|) on 7 in the dark matter halo profile given in 
Eq. We set the other parameters a and /? to be 1 and 3, respectively. The mass of 

the galaxy and the rotational velocity Vh, which are used to determine the parameters 
po and a, are also shown in the figure. 

functions are given by 

dN^ ^ 0.73 e-^-^^- 
dx x^-^ + c/ ' 

where the parameter c/ is 2 x 10~^ for / = W~^W~ and 1.5 x 10""^ for / = 

In the previous studies in Refs. the cutoff parameter c/ is not introduced 
in the fitting functions. However, the behaviors of the simulated fragmentation 
functions at x ^ 10^*^^^'^^ are more moderate than the fitting functions with Cf = 0. 
In this paper we consider cases of heavy EWIMPs (~ 10 TeV) and compare the 
predicted gamma ray fluxes with the EGRET data around 1 ~ 10 GeV. Thus, the 
effect of non- vanishing c/ is not negligible. 

The function N!^r^ for the line gamma ray flux is simply given by 



dK 



(7) 



77 



26{E - m) , (6^ 



dE 

because the dark matter particle is almost at rest in the galactic halo. 
3. Mass density profile of dark matter 

The gamma ray flux depends strongly on the dark matter density profile, because 
it is proportional to the density squared. Many A^-body simulations show that the 
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dark matter halo profiles are given by a universal spherical functional form, 

^^""^ " (r/a)7[l + fr/a)"](/3-7)/" ' ^^^^ 
where a, (3 and 7 are the model parameters. After choosing the model parameters, 
Po and a are uniquely determined by the mass and the rotational speed of the galaxy. 

A famous and frequently used halo model is the Navarro, Frenk and White (NFW) 
profile jH^l, which corresponds to (a,/5,7) = (1,3,1) in Eq. (jU^ . The profile is 
obtained by the numerical A^-body simulation of point particles. Recently higher 
statistical simulations have been performed, and even higher values of 7 are obtained. 
For example, Moore et. al. found the profile corresponding to (a, /?, 7) = (1, 3, 1.4 ~ 
1.5) jHSl- The halo profiles in all A^-body simulations have cuspy structures and the 
density diverges at r = 0. On the other hand, a no-cuspy model has been used for 
a long time, and the model parameters for the King-profile are (a,/?, 7) = (1,3,0). 
It is argued that the rotation curve measurements of low surface brightness galaxies 
disfavor the cuspy profiles [3^, though this disagreement may be resolved by taking 
into account the effect of halo triaxiality [22] • This cusp/core problem is still under 
debate. 

In Fig. ^3 J in Eq. ()66p is shown as a function of the model parameter 7. 
Here we fix a = 1 and (3 = ?>. We use the mass of galaxy interior to 100 kpc 
(Mft(100 kpc) = 6 X lO^^Msun with M^^n the solar mass) and the rotational speed at 
r = 8 kpc (f/i(8 kpc) = 150 km/s) for the determination of po and a in Eq. (|69p. In 
the figure J is sensitive to the value of 7. In the following, we use a moderate value 
J = 500, which is typical for the NFW profile. 

4. Line gamma ray flux from the galactic center 

We discuss the line gamma ray flux from the galactic center due to the EWIMP 
dark matter annihilation by using the cross sections derived in the previous section. 

In Fig. ^[ we show the line gamma fluxes from the galactic center in the cases 
of the triplet and doublet EWIMP dark matters. Here, we take the mass difference 
between the EWIMP and the charged partner as 6m = 0.1, 1, 10 GeV. For the 
doublet EWIMP, 6171^ = 26m is assumed. We also plot the flux obtained from 
the leading-order cross sections in perturbation by broken lines for comparison. In 
Fig. El we show the contour maps of the line gamma ray flux in the MSSM. The 
range of the MSSM parameters is the same as that for Figs. El and [3 in the previous 
section. 
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Figure 12: Line gamma ray flux from the galactic center in cases of the EWIMP dark 
matters. The left (right) figure is for the triplet (doublet) EWIMP. We take the average 
velocity of the dark matter v/c = 10~^, J = 500, which corresponds to the NFW profile, 
and Ail = 10~^. The leading-order cross sections in perturbation are also shown as broken 
lines. 

The large ACT detectors have high sensitivity for TeV-scale gamma rays. For 
example, MAGIC and VERITAS jT3] in the northern hemisphere might reach 
IQ-i^cm-^s-i at the TeV scale while CANGAROO III [HI and HESS [H] in the 
southern hemisphere might reach 10~^^cm~^s~^. From Figs. IT^ and IT^ it is found 
that these ACT detectors may cover broad regions in the parameter space. 

It has been known that the line gamma ray signal is sensitive to the heavier dark 
matter with the mass of the order of TeV, because the annihilation cross section at 
one-loop level is not suppressed by the dark matter mass if the dark matter has the 
SU{2)l charge. Our studies reveal importance of the non-perturbative effects on the 
EWIMP annihilation cross section. After including the effects, the sensitivity of the 
line gamma ray signal to the heavier EWIMP dark matter is enhanced furthermore. 

5. Continuum gamma ray flux from the galactic center 

In addition to the line gamma rays, the EWIMP annihilation produces the con- 
tinuum gamma rays, which come mainly from the decay of in the fragmentation 
of the flnal state particles. Though the number of photons in the continuum gamma 
rays is expected to be higher than that in the line gamma rays, its spectrum may 
lack the distinctive feature. Since the flux of the diffused gamma ray background, 
especially from the galactic center, is not well known, it is difficult to extract the 
annihilation signal (S) from the background (B) as far as S < B. The EGRET has 
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Figure 13: Contour maps of the line gamma ray flux in (M2,/x) planes with tan/3 = 4,40 
in the MSSM. Mi = O.5M2 is assumed in the four top figures while Mi = 3M2 is assumed 
in the four bottom figures. These figures are depicted by using the NFW profile (J = 500) 
and the angular resolution, AO = 10~^. 
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Figure 14: Contour plots of the continuum gamma ray flux from the galactic center in 
the unit of cm~^sec~^GeV~^. The left (right) plot corresponds to the triplet (doublet) 
EWIMP case. Here we take Sm = 0.1 GeV, and other parameters are the same as those 
in Fig. 1121 The shaded regions correspond to (Signal) /(Background) > 1. 

observed the diffused gamma ray emission from the galactic center up to about 10 
GeV [36j. If the power law fall-off of the energy for the diffused gamma ray flux 
^ bg{.E) is assumed, the background B is evaluated from the EGRET result as P, 

"^1^01 = 9.1 X 10-^[cm-^sec-iGeV-i] x ( A(]. (70) 
dE ylGeV / 

In Fig. El we show the contour plots of the continuum gamma ray fluxes from the 
galactic center for the triplet and doublet EWIMPs. Here we fix J = 500, AVt = 10~^ 
and 6m = 0.1 GeV. The shaded regions correspond to 5* > 5, in which B is given by 
Eq. ()70|) . From this figure, it is found that even small regions around the resonances, 
in addition to areas with m ~ 100 GeV, are already constrained from the EGRET 
measurement of the gamma ray flux with the energy 1 ~ 10 GeV. 

In Fig. IT^ regions excluded by the EGRET observation is shown in (M2, /i) planes 
assuming the MSSM. The MSSM parameters are the same as in Fig.^J In this figure 
we use Js with 7 = 0, 0.6, 1 and 1.4. The excluded regions become broader for larger 
7- 

In the near future the GLAST satellite ^U] may observe gamma rays with energies 
in the range 1 GeV^ E^ 300 GeV if the flux is larger than about 10~^°cm~^sec~^. 
The ACT and GLAST detectors will constrain broader regions in the parameter 
space with S > B, or will find the signature of the dark matter. 
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Figure 15: Contour maps of the excluded region by the EGRET measurement of the 
gamma ray flux from the galactic center in (M2,/i) planes in the MSSM for different dark 
matter profiles. Mi = O.5M2 is assumed in the four top figures while it is Mi = 3M2 in 
the four bottom figures, tan (3 = 4, 40. These figures are depicted by using the J factors 
for 7 = 0,0.6, 1, 1.4. We also set the angular resolution to be Ail = 10~^. 
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In this paper we have calculated the pair annihilation cross sections of the 
EWIMP dark matters, which have an SU(2)x charge of the standard model gauge 
group. The leading-order calculation of the cross sections in perturbation is no 
longer valid when the mass is heavy compared to the weak gauge bosons due to the 
threshold singularity coming from the mass degeneracy between the EWIMP and its 
SU{2)l partner (s). The problem have been known for a while for the cases of the 
Wino and Higgsino-like dark matters in the MSSM. We have developed a method to 
take in the singularity and obtain the precise annihilation cross sections. 

We find that if the mass of the EWIMP dark matter is larger than about 1 TeV, 
the attractive Yukawa potentials induced by the weak gauge boson exchanges have 
significant effects on the annihilation processes and the cross sections are enhanced 
by several orders of magnitude due to the zero-energy resonances under the poten- 
tials. As a result, the gamma ray flux from the galactic center due to the EWIMP 
annihilation is enhanced compared to the leading-order calculation in perturbation. 
The line gamma ray flux exceeds the typical sensitivity of the large ACT detectors 
such as CANGAROO III, HESS, VERITAS and MAGIC, IQ-^^^-m) cm~hec-\ in 
the wide range of the MSSM parameters. We also calculated the continuum gamma 
ray flux from the EWIMP dark matter annihilation. The MSSM parameter space 
which is not consistent with the EGRET observation of the gamma rays with 1 GeV 
^ 10 GeV is increased by the non-perturbative effects. The non-perturbative 
effect would also enhance anti-proton and positron fluxes from the dark matter an- 
nihilation. These will be discussed elsewhere |S7] . 

Current observations of TeV-scale gamma rays from the galactic center by the 
CANGAROO and HESS disagree each other in the spectrum. Once it is converged, 
the result may be used to constrain the EWIMP dark matter. 
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